Finite Superposition Solutions for Surface States in a Type of Photonic Superlattices 
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We develop an efficient method to derive a class of surface states in photonic superlattices. In 
a kind of infinite bichromatic superlattices satisfying some specific conditions, we obtain a finite 
portion of their in-gap states, which are superpositions of finite numbers of their unstable Bloch 
waves. By using these unstable in-gap states, we construct exactly several stable surface states near 
various interfaces in photonic superlattices. We analytically explore the parametric dependence of 
these exact surface states. Our analysis provides an exact demonstration for the existence of surface 
states and would be also helpful to understand surface states in other lattice systems. 

PACS numbers: 42.25.Gy, 42.70.Qs, 73.20.At 



I. INTRODUCTION 

Surface states are a kind of localized states found at 
the interfaces between two different media. In 1932, 
Tamm predicted that electronic surface waves may ex- 
ist in a semi-infinite one-dimensional repulsive Kronig- 
Penney (KP) model [l| and these surface states were 
named as Tamm states later. Similar results also exist 
in the attractive KP model @. The energies of Tamm 
states lie in the forbidden gaps for the corresponding in- 
finite KP model. The surface states depend crucially 
on surface termination of periodic potentials It 
has been demonstrated that Tamm states appear when 
the periodic potentials are asymmetrically terminated 
and Shockley states appear while the periodic potentials 
are symmetrically terminated Surface states have 
been observed in several experimental systems such as 
semiconductor superlattices [Q-Q and magneto-photonic 
structures [10| . 

In the past few decades, surface states have been stud- 
ied extensively due to their potential applications in op- 
toelectronic devices. However, an analytical demonstra- 
tion of the existence of exact surface states is still ab- 
sent flTl-[l3j]. For an example, to construct a surface state 
in the semi-infinite one-dimensional KP model, one need 
to calculate analytically or numerically a Bloch wave of 
a complex wave number and then match it with an expo- 
nentially decaying state inside the surface potential \\A\ - 
[l3[ . In a semi- infinite sinusoidal potential |13T | or a semi- 
infinite KP potential [l|, |l4[ , the energies and the wave 
numbers for the surface states are determined implicitly 
by solving a transcendental equation. In addition, there 
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are other approximation methods for constructing sur- 
face states [lj|, such as the coefficient method, the scat- 
tering method, the determinant method and the integral 
method. Up to now, it is still a great challenge to find 
exact solutions for surface states jILC3l — TlSl | . 

Due to the temporal evolution of quantum systems can 
be mapped into the spatial propagation of light waves, 
the engineered photonic lattices provide a highly con- 
trollable platform for exploring similar surface states in 
periodic quantum systems. Surface states near an in- 
terface between a periodic layered medium and a homo- 
geneous medium are found to be analogous to electronic 
surface states in crystals [19( . The formation of Shockley- 
like surface states in an optically induced semi-infinite 
photonic superlattice has been experimentally demon- 
strated [20] • The formation of Tamm states at the bound- 
ary between two periodical dielectric structures has been 
reported [2l[. It has been demonstrated that the nonre- 
ciprocality of the surface modes can be induced by the 
violation of periodicity and the violation of the time re- 
versal symmetry [22j ]. The recent advances on surface 
states in photonic crystals are reviewed in [23j | . 

In this paper, for a kind of infinite bichromatic super- 
lattices satisfying some certain conditions, we find a set 
of the in-gap states in superpositions of finite numbers 
of unstable Bloch waves. These unstable in-gap states 
are then used to construct stable surface states for sev- 
eral typical systems of surface states, such as the semi- 
infinite photonic superlattice, the finite photonic super- 
lattice, two directly coupled photonic superlattices, and 
two indirectly coupled photonic superlattices. The con- 
ditions for the existence of these surface states are de- 
rived. We find that the symmetry and the existence of 
the surface states depend crucially on both lattice param- 
eters and interface parameters. Our analytical results of 
surface states provide an optional benchmark for under- 
standing surface waves in lattice systems. 
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II. IN-GAP STATES IN PHOTONIC 
SUPERLATTICES 



We consider the light propagation in a one-dimensional 
waveguide array. Assuming the waveguide array is 
aligned along the X direction and the light is localized 
along the Y direction, the propagation of the light elec- 
tric field E(X, Z) along the Z direction is described by 
an effective two-dimensional wave equation [24j , 



A dE 
''2^~dZ 



X 2 d 2 E 
'8ir 2 n s dX 2 



U(X)E, 



(1) 



where A is the free-space light wavelength and n s is the 
substrate refractive index. The profile of the effective re- 
fractive index is in form of U (X) = [n 2 — n 2 (X)]/ (2n s ) ~ 
n s — n{X) with the refractive index n(X) for the waveg- 
uide array. By using the periodic modulation techniques 
of n(X) [25|, one can build a bichromatic superlattice of 



f/(X)=n 1 [l-cos(^)]+n 2 [l 



cos(2^ + f 



with the 



amplitudes n\ and 77,2, the modulation period 2 A and the 
relative phase 6. 

By introducing two scaled variables x = irX/A and 
z = irXZ/(4A 2 n s ) and a transformation 4>(x,z) = 
E(x,z)exp[i(n' 1 + n' 2 )z] with n[ = 8A 2 n s ni/X 2 
n' 2 = 8A 2 n s n 2 / X 2 , the system is then described by 



and 
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VW 1 (V, A, x) = exp 
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(exp[— inx] + exp[— i(N — n)x]) 
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Comparing to the Bloch-wave form, the solution 
tj)™ 1 (77, A, x) is the Bloch-wave solution with a complex 
wave number k = N/2 + iA/2. Here, a n and b n can be 
derived from two recursive series: 2r/(n + 2)a n+ 2 + [(n + 
l)(n + 1 — N) + (N 2 + A 2 )/4 - 2 V 2 - ^]a n+l + 2 V {N - 
n)a n + [(n + 1)A - NA/2]b n+ i = and 2rj(n + 2)b n+2 + 
[{n+l){n+l-N) + {N 2 +& 2 )/A-2T 1 2 -f^ 1 ]b n+1 +2i 1 {N- 
n)b n + [NA/2— (n+ l)A]a„ + i = with initial conditions 
of a Q = 1, b = 0, ai = ((3^ + 2 V 2 - (N 2 + A 2 )/4)/2t7 
and bi = — NA/Arj. The propagation constant f3 7 ^ 1 
corresponds to the mi-th real zeros of ^+1(77, A, f3) — 
<iN+i+tt>N+i = in ascending order. Since £^+1(77, A, j3) 
is a polynomial of degree N + 1 in the propagation con- 
stant /3, fi 7 ^ 1 have at most N + 1 solutions. Mathemat- 
ically, one can construct a linearly independent solution 
for Vv 1 in form of t/^ 1 = ip^ 1 ^{ip^^dx. Although 

ip^ 1 an d ^n 1 nave the same propagation constant f} 7 ^ 1 , 
their divergence properties are opposite: ijf^ 1 — > when 
tp^ 1 — ¥ co and vice versa. 



with V{x) — —n' 1 cos{2x) + n' 2 cos(x + 9). Considering 
its stationary states, (f>(x,z) = ijj(x) exp[— i/3z], the am- 
plitude ip[x) obeys a time-independent equation 



(3) 



with P denoting the propagation constant. Obviously, 
the light propagation is equivalent to a quantum particle 
in an external potential. 

By applying the Bloch-Floquet theorem, the solutions 
for Eq. of an infinitely periodic V(x) are Bloch waves, 
ipn,k{x) = exp[ikx)u ni k(x), where k is the wave num- 
ber, n is the band index and u n ^k{x) has the same pe- 
riodicity of V(x). The Bloch waves of real wave num- 
bers are amplitude-bounded oscillatory solutions. Oth- 
erwise, the Bloch waves of complex wave numbers show 
unbounded exponential behavior [29l |30| |. The energy 
spectrum for Eq. ([3]) consists of bands in which there exist 
only amplitude-bounded oscillatory solutions and gaps in 
which there exist unbounded oscillatory solutions. 

Under conditions of 6 = arctan(A/(A r + 1)), n[ = 2r] 2 
and n' 2 = 2rjy/(N + l) 2 + A 2 with integers N > 0, the 
potential V(x) could be denoted by V(r?, A, N, x) and it 
supports a set of in-gap solutions (see more details in the 
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FIG. 1: (Color online). Band-gap structures for V(rj, A, N, x). 
(a) The first band at k — and the in-gap propagation con- 
stant Po f° r N — and 77 = 0.1. (b) The first two bands 
at k — ±0.5 and the in-gap propagation constants fi\' 2 for 
N = 1 and 77 = 0.1. (c) The first two bands and ^ for N = 
and A = 0.3. (d) The first two bands and /3^' 2 for N = 1 and 
A = 0.3. In which, the red dot in (a) corresponds to the one 
in (c), and the two red dots in (b) correspond to the four red 
dots in (d). 



In 



in principle, for a given superlattice V(r],A,N,x), 
one may determine the coefficients for the finite- 
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superposition solutions from the two recursive series. For 
N < 3, one can easily obtain the exact forms for a n and 
b n . However, if TV > 3, it is very difficult to give the ex- 
act forms for a n and b n , and one has to find their values 
by using numerical methods. Below, we consider the two 
simplest cases: N — and N = 1. For the case of N = 0, 
the sole finite-superposition solution is expressed as 



A, a?) = 



exp 



A 

-— x — 2rjcosx 



(5) 



with /3^(r?,A) = -(A 2 + 8r7 2 )/4. For the case of N = 
1, there are two finite-superposition solutions. The first 
finite-superposition solution is in form of 



^(i),A,ie) = exp 
Arj — -y/1677 2 — 
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A fx 
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with /3{(r), A) = (1 - A 2 - 8r/ 2 - 2^/l67] 2 - A 2 )/4. The 
other finite-superposition solution reads as 



Vi(»7, A,x) = exp 



A 

— x — 2?7 cos x 
2 1 



4rj + a/1677 2 - A 2 
477 



/x\ A fx 
cos — — — sin — 
V 2 / 477 V 2 



,(7) 



with /3 2 (?7, A) = (1 - A 2 - 8?7 2 + 2 V / I6r ] 2 - A 2 )/4. If 
I677 2 — A 2 =0, the two in-gap waves ipi(r), A, x) and 
V? (77, A, x) are identical. 

In comparison with the band-gap structure, if A / 0, 
we find that (3q falls into the semi-infinite gap below the 
lowest band and (3{ lies in the first band-gap, see Fig.Q] 
This means that these finite-superposition solutions are 
a kind of in-gap states. If A = 0, the finite-superposition 
solutions become stable Bloch-wave solutions, since the 
wave numbers become real. Interestingly, P\ and /3 2 form 
a closed circle connecting the first two bands at A = 0, 
see Fig. \T[b). Moreover, the in-gap state Vo appears 
at k = 0, while the in-gap states Vi' 2 appear at k — 
±0.5, see Fig.QTc) and (d). Since the in-gap states grow 
without bound, they are unphysical states for the infinite 
periodic system. However, as we will show below, the in- 
gap states can be used to construct a special class of exact 
surface states in several typical models. 



III. SURFACE STATES IN SINGLE-INTERFACE 
SYSTEMS 

One of the most famous single-interface systems is a 
semi-infinite periodic system of a truncated V(v, A , N, x) 
connecting a constant refractive index Vq |ll| - [l4| . The 
potential for such a system reads as, 



Vo, x < x (region I), 

V(r), A, N, x), x > xq (region II). 



For an allowed surface state, in the region II, it should 

be in form of V™ 1 (v, A, x) for A > (or V^ 1 for A < 0), 

Vn(x) = C 2 V7 (or Vn(z) = C 2 V>)- Below we will 
only consider the case of A > 0. In the region I, ipi(x) = 



P^x] if Vo > Pn 1 - The coefficients d are 



d exp[\/Vb - HN 
determined by the normalization condition. By applying 
the continuity condition at the interface x = xq, we find 
that 



^V Q -(3^( V ,A)^Wp(r h A,x ), 



(8) 



with W^( V ,A,x) = Vr(^A,x)/V7(»?,A,x). Here, 
the dot denotes the derivative with respect to x. Thus 
the surface state exists if W'^ 1 (77, A, N, xq) > and 
Vo = P^ 1 + (Wp) 2 . For example, in the simplest 
case of N = 0, the interface parameters xo and Vo sat- 
isfy the conditions, sin(x ) > A/4?7 and Vo — —(A 2 + 
8t7 2 )/4+ (A/2 - 277sin(x )) 2 . In Fig.^a), we show the 
surface wave in this semi-infinite periodic system with 
7/ = 0.3, A = 0.2 and N = 0, which corresponds to n-y = 
2.2 x 10 -4 , 7i 2 = 7.6 x 10~ 4 ) and 6 = arctan(0.2) in an ex- 
perimental system of A = 980 nm, n s = 1.518 and A = 8 
/jm @,[2^|. In this situation, we take xo = n/2 and thus 
have V = -(A 2 + 8t7 2 )/4 + (A/2 - 277sin(x )) 2 = 0.06. 




FIG. 2: (Color online). Surface states in single-interface sys- 
tems, (a) Vi(?7, A, N, x) of x = tt/2, r\ = 0.3, A = 0.2, 
Vo = 0.06 and N = 0. (b) V2 (t?, A, N, x) of x = tt/2, 77 = 0.1, 
A = 477 sin xo and N — 0. 

Another typical single-interface system is of two trun- 
cated periodic potentials connecting at the interface. We 
consider a system of V(— 77, — A, N, x) and V(rj, A, N, x) 
with A > connecting at the interface x = xq hrj. its 
potential is expressed as 



V 2 (r,,A,N,x) 



V(-T], -A, N, x), x < x (region I), 
V(r], A, N, x), x > xq (region II). 



Thus, we have Vi(x) = CiVat 1 (— f], — A, x) for x < xo 
and Vn(x) = C2V/V 1 (.Vi A, x) for x > xo- The continuity 
conditions at x = xo give 

(-77, -A, N, xo) = (rj, A, N, x ). (9) 
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For the case of N = 0, it is easy to find 4 77 sin xq = A from 
the continuity condition at the interface. Therefore, the 
surface wave exists only when |A/4?/| < 1. In Fig. EJb), 
we show the surface state for N = and xq = ir/2. 
We also find that the position of the interface xo affects 
strongly the shape of the surface waves. For an example, 
the surface wave for the case of xq — 7r/2 is symmetric 
about x = xo = 7r/2. While the surface wave for the case 
of xq = 7r/6 is asymmetric about x = xq = 7r/6. 



IV. SURFACE STATES IN 
DOUBLE-INTERFACE SYSTEMS 

One typical double-interface system is a finite periodic 
system V(r), A, N, x) sandwiched by two constant refrac- 
tive indices Vq and V\ [3||, which obeys the potential 

Vo, x < x (region I), 

V 3 (rj, A, N, x) = { V(r], A, N, x), x < x < x x (region II), 
Vi, x > x\ (region III). 

In the region II, one may use the finite-superposition 
solution ipN 1 (Vi A, x ) and its linearly independent solu- 
tion ip 7 ^ 1 to construct the surface state, that is, ipii{ x ) = 
C2VV 1 ("Hi A, x) + C2ip^ 1 ■ In other two regions, the physi- 
cal state must be non-divergent and normalizable. There- 
fore, if Vb,i > P™ 1 , we have ip^x) = d exp[y/V - ffix] 
for x < x and ipm(x) = C 3 exp[— ^/Vi — j3^x] for 
x > x\. Similarly, the continuity conditions at the two 
interfaces x — xq and x = x% request 



mi 
N 



R 



1 + R 



(10) 



with 



R — C2I C2 — 



y/V -ffi-Wp(T,,A,X ) 



ami prai 

Pn r 



N 



Oo)' 



KT(x) 



and 



FT(x) 



In principle, the ratio between C2 
and C2 can be arbitrary. However, to satisfy the 
continuity conditions at the two interfaces, the interface 
parameters (Vo, V\, xq, x\) and the lattice parameters 
(77, A, N) should obey some certain conditions. In Fig. [3l 
we show two surface states for 77 = 0.3, A = 0.1, N = 
and xo = —x\ = — l9n/2. If R = 0, the continuity 
conditions request Vq = 0.12 and Vi = 0.24 and the 
surface state is localized around the interface x = xo , see 
Fig. (3Ja). If R = 00(6*2 = 0), the continuity conditions 
request V = 0.196151 and V x = 0.142676 and the 
surface state is localized around the interface x = x\, 
see Fig. [31(b). If R = 0.0696, the continuity conditions 
request V = 0.123884 and Vi = 0.147421 and the 
surface state is almost equally localized around x — xq 
and x = Xi, see Fig. [21(c). 
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FIG. 3: (Color online). Surface states in a finite periodic 
system described by VJs(?7, A, N, x) with N = 0, 77 = 0.3, 
A = 0.1, and x = -xi = -19tt/2. (a) Vo=0.12, Vj = 0.24, 
and R = 0. (b) Vo = 0.196151 and Vi = 0.142676, and 
R = 00. (c) Vo = 0.123884, Vi = 0.147421, and R = 0.0696. 



Another typical double-interface system is a constant 
refractive index Vq sandwiched by two truncated periodic 
systems V(rj, -A, N, x) and V(r], A, N, x) ^MM- The 
corresponding refractive index profile is in form of 

V(r), -A, N, x), x < x Q (region I), 
Vi(?7, A, N, x) = ^ Vo, x < x < Xi (region II), 

V(r}, A,N,x), x > xi (region III). 



Thus, we have ipi{x) = (?7, — A, x) in the 

region I, ipn(x) exp[- A /Vb - /3" r ll x 

C7 2 + exp[+ x /V r ; 



i^m(x) = Ca^^r?, A, x) in the 
continuity conditions request 



C 2 exp[-^Vb - , 
in the region II, 
region III. 



and 
The 



W^faA,^) _ R-exp[2y/V -^ Xl }-\ 



ami 
Pi* 



R-exp[2 y /V -l3S 1 x 1 ]+r 



(11) 



it.h 



R 
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^Vo-PT +W% 1 (y,-A,x ) 
vW-^(lrA,io) 
xexp[-2 x /Vo-/3™ 1 x ]. 



In Fig. 01 we show two surface states for 77 = 0.3, A = 0.1 
and N — 0. The surface waves strongly depend on the 
two interface positions. For xq — — 7r/2 and x\ — ir [or 
x = — 77r/6 and x\ — 7r/2], the surface state has an 
asymmetric distribution, see Fig. Ufa) [or (b)]. While for 
xo = — 7r/2 and x\ = 7r/2, it becomes symmetric, see 
Fig. He). 



V. CONCLUSION 

In conclusion, by using the superpositions of finite 
numbers of unstable Bloch states for the correspond- 
ing infinite periodic systems, we have given an efficient 
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FIG. 4: (Color online). Surface states in a double- interface 
system described by V4(?7, A, N, x) with r\ = 0.3, A = 0.1 
and N — 0. (a) xq = —n/2, xi — n, Vo — 0.13, and R = 
0.0255412. (b) x = -7tt/6, x x = tt/2, V = 0.120845, and 
R = 254.28. (c) x = -xi = -tt/2, P = 0.30, and R = 1.0. 

approach for constructing a special class of exact sur- 
face states in photonic superlattices. These exact surface 
states have the same propagation constants for the finite- 
superposition states in the energy gaps and so that they 



are a kind of stable in-gap states. This method has been 
used to find parts of the exact surface states in several 
typical systems involving a single or two interfaces. By 
matching two solutions (a finite-superposition solution 
and a free-space solution or two finite-superposition so- 
lutions) at two sides of the interfaces, the existence con- 
ditions for surface states are obtained analytically from 
the continuity conditions. The existence and the shapes 
of the exact surface states not only depend on the inter- 
face parameters, but also rely on the lattice parameters. 
Our results give an analytical demonstration of the ex- 
istence of the surface states and should shine light on 
understanding and controlling the surface waves. 
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Appendix: Derivation of the Bloch-wave solutions 

Vv 1 (v, A , x ) 

Below, we give more details about how to derive the 
in-gap solutions 



fpN 1 ^^^) = ex P 



N A . 
i \ — + i— \ x — 2r] cos x 



N 

E 



niPfj ) ( ex pj_j na .] _|_ exp[— i(N — n)x]) 



n=0 



(exp[— inx] — exp[— i(N — n)x]) 



(A.1) 



for the Schrodinger equation 



dx 



— i>{x) + V(x)1>{x) = W(x), (A.2) 



with 

V(x) 



-n'i cos(2a;) + n 2 cos(x + 9) 

-n'i cos(2a;) + n' 2 cos 9 cos x — n 2 sin 9 sin a;. 



In general, the in-gap solutions stay at the edges of 
Brillouin zones and the imaginary parts of their wave 
numbers are continuous in a finite region [29L l30l |. To 
find the explicit expression for a specific set of in-gap 
solutions, we apply the following transformation 



with 



ip(x) 



A = 



exp[— ix\ 



exp[ 



2rai cos x]f A </>(£)' 



(A.3) 
(A.4) 



' 2n'i — n' 2 cos 9 — in' 2 sin 9 



2 v /2^[ 



Apparently, the parameter —A denotes a complex wave 
number. From the Schrodinger equation (A.2), we have 



,d 2 4> 



2n[ - V2nir + (2A + 1)£ 



d<p 
d£ 



A 2 - /3 - ni - (y^n{-n 2 cos ej^ 



0. 

(A.5) 



By writing the solution for </>(£) as a standard power- 
series expansion, 



(A.6) 



with the initial condition do = 1 and d-i = 0, where 



Cq(u) — ri 2 cos 9 — yj2n' l {n + 1), 
ci(n) = (n + l)(n + 2A + l) + A 2 
c 2 (n) = y/2n[(n + 2). 



Usually, the solution (A.6) is an infinite series. How- 
ever, similar to the procedure of obtaining the Hermite 
polynomials for a harmonic oscillator, one can impose the 
truncation condition dj — with j > N + 1 and then the 
solution (A.6) becomes a finite series. If d/v+i = and 
An +2 = 0, we can get all following dj = with j > N + 2. 

To obtain d^+i = 0, from the three-term recurrence 
relation (|A.7P with n = N — 1, we have 



c (N-l)d N -i + Cl (N-l)d 



N 



0. 



(A.8) 



For a given periodic lattice, due to <i/v is a polynomial 
of degree N in /3, this equation requests that the propa- 
gation constant /? must be a solution for a polynomial of 
degree N + 1 in fi. 

To obtain d/v+2 = 0, from the three-term recurrence 
relation (IA.7I) with n = N, the coefficient cq(N) should 
satisfy 



co (AO 



id - \f2rif x {N + 1) = 0. 



(A.9) 



Clearly, this equation requires a special relation between 
the lattice parameters n' ly n\ and 9. 

Therefore, under the conditions (IA.8I) and (|A.9[) . the 
series solution </>(£) becomes a polynomial, 



<Me) 



N 

E 

71 = 



(A.I0) 



After some mathematical calculation, we get the follow- 
ing in-gap solution 



ip{x) 



exp 



N . n' 2 sin 9 \ 
~2 +l 2yf2n T 1 J X 



'2n' x cos x 



N 



x^>„(/3 
n=0 



) exp I 



(All) 



Here, j5 1 ^ 1 denotes mi-th real zero of d/v+i = in ascend- 
ing order. It is clear that ip(x) are Bloch-wave solutions 
with the complex wave numbers k = —A = N/2+i n z sln ® , 



If ^/2n\ = 2-q, n' 2 cos 9 = 2n{N + 1) and n 2 sm9 = 2nA, 
we have n[ = 2n 2 , ri 2 = 2n^(N + l) 2 + A 2 and 9 = 
arctan(A/(iV + 1)). The complex wave number is given 



one can easily find that the coefficients d n are determined 
by a three-term recurrence relation, 



k = 



N 



.A 



(A.12) 



Co{n)d n + ci{n)d n+1 + c 2 {n)d n+2 = 0, 



(A.7) 



Due to V(x) is a real function, we can take the real 
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part of i/j(x) as a solution of Eq. (|A.2[) 



N A . 
i | — + z — ) x — 27/ cos x 



N 

E 



a "(ftv ) ( ex p[_j na; ] _|_ exp[— i(iV — n)x]) 



bJj3 mi ) 

+i— — — — (exp[— inx] — exp[— i(N — n)x\) 



Given c (0) = 2r?,ci(0) = -(A + i) 2 /4 - 2^ 2 - /3 and 
rfi = ((3 + (A - i) 2 /4 + 2T7 2 )/2r7, we have 

(jfl + (A + z) 2 /4 + 2? ? 2 )(/3 + (A - «) 2 /4 + 27? 2 ) - 4t? 2 = 0. 

(A.17) 

This equation have two real roots, Pl(f], A) = (1 — A 2 — 
8r, 2 - 2^/l6?? 2 - A 2 )/4 and /3?(r?, A) = (1 - A 2 - 8t? 2 + 
2y/l6rj 2 - A 2 )/4. For ^(77, A), we have a = 1, b = 
Q,ai = — -\/l6r/ 2 — A 2 /4?7, and b\ = —A/4i], therefore 
(A. 13) the corresponding in-gap state reads as 



where c^/?™ 1 ) = M/3™ 1 ) +ifr„(/3™ 1 ). This completes 
the derivation of Bloch-wave solutions (|A.1[) . 

In the following, we show how to give the in-gap waves 
for the cases of N = and N = 1. For the case of N = 0, 
from Eq. (|A.8|) with AT = 0, we have 



ci(— l)do = A 2 — rij — 
A 2 

= -— - 2r? 2 - = 0. 



(A. 14) 



This equation has only one real root (3q = — (A 2 + 8?7 2 )/4. 
From Eq. (lA.lj) with N = 0, the corresponding in-gap 
state is given as 



ipoiv, A,x) = exp 



-x — 2?7Cosa; 



(A.15) 



For the case of AT = 1, from Eq. (|A.8[) with iV = 1, we 
have 



ipifa A,ar) = ex P 



A 

x — 2ri cos x 

2 ' 



4t; - Jl&if - A 2 / x \ A 

- cos — — — sin 

4rj V 2 / 477 



(!) 



(A.18) 



Similarly, for j3 2 {rj, A), we have ao = 1, 6q = 0,ai = 
-\/l6?7 2 — A 2 /4?7, and bi = —A/477. The corresponding 
in-gap state reads as 



i>i(v, A, 1 ) = CX P 



A 

— x — 2r7C0sa; 
2 ' 



477 + yjl&ri 2 - A 2 
477 



(!)-£-(!) 



(A.19) 



c o (0)d o + ci(0)di =0. 



(A.16) 



